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ABSTRACT 

We study the effects of radial flows on Galactic chemical evolution. A simple analytic 
i-C ■ scheme is developed prescribing the coupling of infall from the intergalactic medium 

and radial flows within the disc based on angular momentum conservation. We show 
that model parameters are tightly constrained by the observed [Fe/H]-abundance gra- 
dient in the Galactic disc. By this comparison the average rotational velocity of the 
onfalling material can be constrained to 0.7 ^ v/V c ^ 0.75, or respectively ~ 160km/s 
when assuming a constant disc circular velocity of V c = 220km/s. We test the robust- 
ness of this value against the influence of other processes. For a very simple model 
of inside-out formation this value changes only by Av/V c ~ 0.1, i.e. ~ 20km/s, and 
significantly less on more realistic scenarios, showing that inside-out formation does 
not alone explain the abundance gradient. Effects of other uncertain parameters, e.g. 
star formation history and star formation efficiency have very small impact. 

Other drivers of inflow beyond our explicit modelling are assessed by adding a 
fixed inflow across the whole disc. The churning amplitude only mildly affects the 
results mostly by slightly flattening the metallicity gradient in the inner disc. A new 
process causing radial gas flows due to the ejection of material by stars moving on 
non-circular orbits is studied and seems to contribute negligibly to the total flows. 

We further show that gaseous outer discs cannot be the main source feeding the 
persistent star formation in the inner regions by a direct inflow. 

Key words: - Galaxy: evolution - Galaxy: abundances - (galaxies:) intergalactic 
medium - galaxies: ISM - galaxies: kinematics and dynamics 



1 INTRODUCTION 

Radial flows of gas have been studied extensively in mod- 
els of chemic al evolution. Thei r impo rtance has first been 
recognised bv lTinslev fc Larson! [|l978u . who emphasise their 
effects on galactic abundance gradients and conclude that 
these cannot be ignored in any viable model. lLacev fc Falll 
1 19851 ) discussed three main drivers of these flows. First, on- 
falling gas might have a different specific angular momentum 
than the gas rotating in the disc. As the onfalling gas is sus- 
pected to have a lower rotational velocity, these flows are 
expected to be directed inwards at a velocity in the range of 
a few km/s l|Lacev fc Falllll985l ). Second, friction within the 
rotating gas layers might cause radial flows. A viscous disc 
is subject to an angular momentum flux whenever the angu- 



lar speed changes with radius. A classic Galactic disc with 
its nearly flat rotation curve has a global outwards directed 
angular momentum flux, which drives inflows in the inner 
regions and an expansion of its outskirts. Est imates of this 
effect yield velocities of the order of 0.1 km/s (|Lacev fc Falll 
1985). Third, interactions of the gas with a central bar and 
spiral patterns, especially around the corotation resonance, 
will cause radial flows, t ypically directed awa y from the po- 
sition of the resonance. lLacev fc FalU l|l985h cite velocities 
of 0.3 km/s at a distance of lOkpc from the galactic cen- 
tre. Further the gas loses angular momentum to the stellar 
population by interactions/shocks along the spiral arms. A 
comparison of these quantities shows that especially stronger 
flows are preferentially driven by angular momentum dilu- 
tion via onfall. 
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Infall has been kno wn to affect t he evolution of disc 
galaxies for many years (|Larsorj|l972h . Even if no gas at 
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all was lost to the intergalactic medium, roughly one third 
of the mass involved in each star formation cycle would get 
locked up in stellar remnants or very long-lived low mass 
stars and the star formation in spiral galaxies would ex- 
haust the available gas reservo irs on cosmologically short 
timescales l)Larson et al.l Il980f ). Hence the observed star 
formation rates require ample infall of fresh gas. Infall is 
also required from a chemical evolution point of view to 
avoid exceedingly high metallicities after a few Gyrs and 
to avoid deuterium abundances far below the observations 
|Linskv et al.ll2006l) . The mode of this accretion is obscure: 
It is frequently adopted that gas is delivered directly to 
alaxies by infalling satellites. This direct infall is observed 
Sancisi et al.ll2008l ). but they also find that it brings about 
an order of magnitude too little material to sustain the cur- 
rent star formation rates. Further, such "cos mological infall" 
leads to rather unsatisfacory outcomes (see IColavitti et al.l 
2008) in chemical evolution models, since the clumpy accre- 
tion and star formation produce unobserved excursions in 
the detailed abundances, especially in the [Fe/H] — [a/Fe] 
plane. 

Analytical chemical evolution models assumed a vanish- 
ing a ngular velocity of onfalling material (|Mavor fc Vigrouxl 
Il98ll ). Later onfalling material with a constant fraction of 
the disc's angul ar velocity at the g i ven r adius was con- 
sidered as well (|Pitts fc Tavlerl 1 19891 1 19961) . But no argu- 
ment aside from trying to fit the observed abundance gra- 
dients could be made for any particular prescription. In a 
different line of development radial gas flows were imple- 
mented to classical chemi cal evol ution models, starti n g with 
iGoetz fc Koeppenl J 19921). later iPortinari fc Chiosil (|2000l ) 



and recently bv lSpitoni fc Matteuccil (|201ll ). Those models 
assume preset patterns of the inflow velocity in the disc, but 
do not draw a significant link to physical properties like the 
angular momentum budget. 

It seems natural, though, to link the available angular 
momentum budget to the puzzle of Galactic accretion. The 
deficit of visible cold gas infall was recently connected to 
the presence of a concurrent mode of delivery: Hot coro- 
nal gas can be swept up by turbulent cooler gas clouds 
that are driven out of the Galact i c plane by supernova 
feedback (Fratcrnali fc Binnevl I20061 . |200S| ; iMarinacci et~al1 
|2010| ). Dynamic arguments and observations strongly sug- 
gest that the coro nal gas should lag be hind the rotation 
speed of the disc (jMarinacci et al.ll201ll ). Theoretical cal- 
culations for the Milky Way predict values for the global 
accretion rate sufficient to sustain its star formation and 
put an estimate of the lag of accreted coronal material at 
around 75 km/s compared to th e disc rotational velocity of 
220km/s |Marasco et aljfeoill ). Hence in contrast to cos- 
mological cloud accretion this offers a firm framework with 
relatively well constrained properties of accretion. In addi- 
tion they derive a specific profile for the accretion of gas 
onto the Galactic disc. 

A process neglected so far is the systematic lag of stel- 
lar yields themselves behind the local gas velocities as well 
as friction between the stars and the interstellar medium. 
This is mediated by the larger velocity dispersion and hence 
asymmetric drift of stars. We will discuss this in the follow- 
ing paper. 

The purpose of this paper is to give a better parametri- 
sation of flows than that used originally in the model of 



SB09 in the sense that it provides sufficient freedom to fit 
the data and depends on physically meaningful parameters. 
The model crucially relies on the assumption of conservation 
of angular momentum which may not be true in all cases, 
but which should be a plausible first approximation to un- 
derstand the overall flows in the Galaxy. Our approach is 
different from those found in the literature as we assume a 
radially varying profile for the angular velocity of the on- 
falling gas. Fitting the model to the data provides informa- 
tion on the kinematic properties of the onfalling gas. These 
properties depend on the proposed origin of the onfalling 
material and the assumed way of accretion onto the Galaxy. 
A comparison of the predictions derived from a model of 
chemical evolution with those of theories for the origin of 
the onfalling material might help to shed light on the ques- 
tion of how galaxies accrete gas. 

The paper is organised as fol lows. Section[2]summarises 
the aspects of the SB09 model (|Schonrich fc Binnevl \200$ ) 
important to this work. In section[3]we give the equations of 
the new flow model to be used in the simulations. We model 
the angular velocity of the onfalling gas as a linear function 
of galactocentric radius. The general behaviour of the model 
and the dependence on its parameters is presented in sec- 
tion [3] In section we fit the new model to the Galactic 
radial abundance gradient to constrain the physical model 
parameters. We inve stigate the imp li cation s of the accretion 
profile proposed by iMarasco et al.l (|201ll ) in section \E\ In 
section [7] we study the proposed new mechanism for caus- 
ing radial flows due to the lag of stellar yields. Section [8] 
discusses an implementation of inside-out formation in the 
model and section [9] explores limits on the accretion from 
an extended gas disc. Section \W\ sums up. 



2 THE UNDERLYING MODEL 



As a s tarting point we use the model oflSchonric h fc Binnevl 
(2009) (SB09). In the following section we briefly recapitu- 
late those parts of the model that are of particular impor- 
tance to this work. 



2.1 The setup 

The simulation follows the evolution of a disc of size 20 kpc 
divided into 80 annuli, which are equispaced with a width 
of 0.25 kpc. The model has discrete timesteps of 30 Myrs. 
Each annulus contains a cold gas phase, which produces 
stars, and a hot phase, which is fed by stellar yields and 
condenses back onto the cold phase. The model tracks the 
detailed abundances in each element considered and holds 
information over the stellar populations formed at each 
timestep. 

2.2 Star formation rate and initial mass function 

The star form ation rate (SFR ) is assumed to follow a simple 
Schmidt law |Schmidtlll959l ). i.e. it depends on the surface 
density of the gas (E s ) according t o a p ower law, whose 
index was determined bv lKennicu tt (1998). We hence use 



dt 



C else, 
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where C is chosen such that the SFR is continuous in E and 
k is fitted to observations. For simplicity we assume that 
the critical density E c = throughout this work. A positive 
value would produce a star formation break in the Milky 
Way, beyond which the star formation rate drops sharply. 
Such a break should not take place out to 12 kpc and be- 
yond this point it has only very minor effects on the local 
chemistry. 

For the initial mass function (IMF) the model uses a Salpeter 
function with upper and lower mass limits of 0.1 Mq and 
100 M Q 

dN , ,_ 2 35 

— oc M 235 . 
dM 

Apart from some influence on detailed abundances (which 
are of no interest for this work) by the shape of the IMF at 
high solar masses, it mainly impacts the lock-up of mate- 
rial at lower masses and by this the need for fresh material 
to sustain a certain star formation rate and influences the 
general metallicity level in the disc. In good agreement with 
general findings we assume that there is at best a weak de- 
pendence of the IMF on Galactocentric radius. 



2.3 Flows and onfall 

The disc is initialised in the first time-step with a gas mass 
of Mo and an exponential surface density profile. Subsequent 
onfall increases the gas mass of the disc by 



M 



Mi 



-t/bt 



, M 2 . 



t/b 2 



Starting from a relatively low initial mass Mo, the double 
exponential law allows for a rapid build-up of the Galactic 
disc mass and SFR on the short timescale bi, followed by a 
slower decline set by the long timescale 62. Such an onfall 
law with two timesc ales has been widely used in Galactic 
models starting with IChiosil (|1980T ). In contrast to models 
without radial flows, in which the radial structure of onfall 
can be inferred from the present mass densities, only the 
total accreted mass can be constrained in our model at this 
point. 

In general the mass balance at each ring then reads 



Mi = (Oi - Di) + (Ti - *i) + A, 



(1) 



with onflow O, outflow/loss to the IGM or to the corona 
D, star formation 9, stellar yields T and the redistribution 
balance A from mass exchange with neighbouring annuli. 

To preserve an exponential disc profile we hence have 
two free functions to fulfill a single constraint: The spatial 
distribution of infall and the redistribution of material 
throughout the disc. To deal with this underdetermination 
SB09 used a very simple parametrisation consisting of the 
two parameters f a and fb, which determine the shape of 
radial inflow and infall from the IGM. As a main ingredient 
of this paper we will replace this simple parametrisation 
with the far more physical condition of angular momentum 
conservation and thus trace the redistribution of material 
back to angular momentum exchange and the specific 
angular momentum of onfall. 



2.4 Kinematics 

To cope with Galactic dynamics and kinematics, the model 
includes a machinery to calculate the spatial distribu- 
tion and kinematic of stars using a radius-dependent age- 
dispersion relationship under the assumption of isothermal 
subpopulations. More importantly it parametrised radial mi- 
gration in gas and stars via an effective exchange of mass 
b etween different annuli in co ncordance with the mechanism 
of lSellwood k, Binnevl (|2002T ) and by reproduction of the lo- 
cal metallicity distribution proved the existence of signif- 
icant radial migration in our Galaxy. The mass exchange 
happens with adjacent rings at each sub-timestep, allowing 
for two exchanges per time-step. The parametrisation de- 
pends on the churning amplitude, a free parameter of the 
model, which in this work we keep at the preferred value 
of k c h = 0.25 if not explicitly stated otherwise. In any case 
slightly stronger or weaker mixing of the gas and stellar 
yields has a minor impact on the metallicity gradient of the 
gas. The value of SBI was derived on the assumption a ra- 
dial abundance gradient that matched a blend of different 
datasets throughout the Galactic disc. This gradient domi- 
nates the width of the local abundance distribution, so that 
a shallower gradient needs to be balanced by higher churning 
rates. 



3 NEW FLOW MODEL EQUATIONS 

In this section we outline how we drive radial flows through 
the Galactic disc by angular momentum conservation with 
respect to the onfall of fresh material from the intergalactic 
medium (IGM) or the Galactic corona. 

In the following we will assume a flat rotation curve, i.e. 
a fixed circul ar speed V c . While this is anyway a good ap- 
proximation |Foster fc Cooper] I2010T I. moderate deviations 
will have very weak impact on our results, espcially since 
we will define the velocity of infalling material in units of 
the local circular speed. We further assume in concordance 
with the observations that the gas is on nearly circular or- 
bits. 

Three functions in Galactocentric radius R are of im- 
portance to the infall problem: The surface density of gas 
already present in the Galactic disc E p (i?), the specific an- 
gular momentum or respectively azimuthal velocity of the 
onflowing/accreted material v (R) and the surface density 
flux of the onfalling material s (R). E P (_R) is essentially 
known from the model history, while there are hardly any 
constraints on v (R) and s (R). Only the total amount of 
infall is preset by the model and hence the problem is un- 
derconstrained by the only constraint of the surface den- 
sity profile after infall. As in SBI we can now use measured 
abundance gradients to constrain the functions. To keep the 
number of free parameters as low as possible while allowing 
for some freedom we assume a linear dependence of v (R) 
on the Galactocenric radius and parametrise 



v„(R) = V c (b + a-^- 

\ Rout 



(2) 



where we assume an outer disc radius R ou t = 20 kpc fit to 
the limit of the calculation, V c is the circular speed of the 
Galaxy, b is the central (R = 0) relative angular momentum 
of the accreted material and a + b is the relative angular 
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momentum at the outer edge. Once v (R) is known, the 
onfall density profile s (R) is fully determined by the aimed 
at gas surface density. 

In the following we will typically consider the case that 
the angular momentum of onfalling gas is lower than that 
of the circular orbit and speak of radial inflow even if radial 
outflow would occur if the angular momentum of the on- 
falling material were higher than that of the disc gas. This is 
not strictly necessary for the general formalism introduced 
here, but covers the parameters relevant to this work and 
allows a simple scheme of solving the equations iteratively. 

For convenience we integrate the onfall over a small time 
(our timestep) to obtain S (iZ) the surface density of mate- 
rial dumped over that time onto the disc. In each annulus 
this dilutes the angular momentum causing a deficit in rel- 
ative angular momentum AJ compared to a circular orbit: 



AJ a = 2tt 



(v (R) - V c (R))RT, a (R,t)RdR, 



(3) 



where Vd(R) is the rotation curve of the Galaxy and T, (R, t) 
is the area density of the onfalling gas. Assuming that the 
surface density of onfalling material is a smooth function 
and considering the relatively fine spacing of the model, we 
can well approximate the onfalling mass onto each ring as 
M 0yn ~ 7rS (ii^ +1 + i? 2 ) by neglecting the spatial variation 
of E over the bin. Moving this out of the integral, we hence 
obtain: 



= 2 M ,n(fi,t) f R ^ {vo(R) _ MR)) R l dR (4) 
rt n+l -"71 J R„ 



With the assumed flat rotation curve Vd(R) = V c and using 
the parametrisation of equation ([2]) the integral solves to: 



/' 

J R, 



(v (R)-V c )R 2 dR- 



i — Rn) + 



'' 1 (R 3 n+1 -R 3 n )) (5) 



Setting an equidistant spacing with width AR for the N 
rings, we get 



AJ = 



2a 4n 3 + 6n 2 + 4n + 1 



+ 



4iV 2n + 1 

+ 2(6 - 1) 3n 2 + 3n + 1 



3 2n + l .-^ »""<■»■ ! <« 

Now we can make up the angular momentum balance com- 
pensating the angular momentum loss caused by onfall via 
an inflow through the disc 

= AJ tot (R) = AJ {R) + AJ l (R), (7) 
where 

AJ,(R) = Mi(R,t)V c AR (8) 

is the angular momentum given up by the radial flow of disc 
gas and Mi (R, t) is the mass flowing in from the next ring. 
Plugging equation ([6]) into equation (0 yields 



Mi(R,t) 



M a (R,t) \4N 



-2a An d + 6n 2 + An + 1 



2n 



+ 



2(1 - b) 3n 2 + 3n + 1 
3 2n + 1 



(9) 



One might be concerned about the explicit factors of AR and 
N possibly making the results of the simulation resolution 
dependent. However, for our choice of the angular momen- 
tum of onfalling g as equation (|9l) is eq uivalent to expression 
7 in the paper bv lLacev fc Fall \ 19851 ) as can be easily seen 
by using the relation between masses (M /i, respectively for 
onfalling and inflowing gas), their surface densities (a , a g 
for onfalling and disc gas) and the inflow velocity (y>i). 



Mi(R,t) = 2na g RviAt 
M a (R,t) = 2ttE RAR 

Vi(R,t)= 



<7 R(v - V d ) 



2-KOgRAt 



I'd 



(10) 
(11) 
(12) 



Holding a discretised version of these expressions, we test 
that our result is not affected significantly by resolution 
problems in the appendix, section [TT] 

Equation (|9} couples flows and infall by the conserva- 
tion of angular momentum. This coupling combined with 
the total need for fresh gas in each ring determines the total 
amount of mass falling onto each ring from the IGM and 
radially flowing from the next outer ring. Here we consider 
exclusively the case of lower angular momentum of the on- 
falling gas, i.e. the gas will always flow radially inwards. As 
in this scheme each ring just receives mass from its outer 
neighbour, it can be solved in a simple way running against 
the flow direction, i.e. iteratively from the centre to the out- 
ermost ring (like the scheme of SBI). This procedure is only 
viable as long as the flow has a single direction throughout 
the whole Galaxy, otherwise a more sophisticated method 
would be required. The outermost ring is assumed to be 
supplied by a gas layer extending beyond that of the stellar 
disc modelled with the metallicity of the IGM. 

The presented scheme fails if the mass of the next outer 
ring is not sufficient to cover the required inflow. In such a 
case we limit the inflow to the total mass of this ring, which 
is equivalent to increasing the momentum of the infalling 
gas. This procedure is somewhat arbitrary, and the prob- 
lem could be straight-forwardly tackled by finer timesteps. 
Yet, it sets a natural limit to the maximum flow velocities 
of ~ 8.2km/s within the constr aints of 5-10km/s from ob- 
servations of external galaxies (|Wong et al.ll2004h . Besides 
being a plausible limit, at our examined parameter space it 
bears relevance only at the earliest timesteps of the simula- 
tion when the disc builds up rapidly. For the inflow model 
to make sense, this adjustment should not be necessary for 
later times in the simulation. We made sure that our choice 
does not affect the results as compared to a finer time reso- 
lution. 

In the common case onfalling gas should in all annuli 
have less or equal angular momentum than the disc gas, 
which constrains a + b ^ 1. Of special interest are the case 
of constant azimuthal velocity a = 0, which has been studied 
in the literature, and the case of rigid rotation with b — 0. 
We also point out that a fixed relative angular momentum 
of onfalling material is by no means synonymous to a con- 
stant inflow speed in the disc: Even if the onfalling mass has 
the same radial density distribution the same inflow speed 
needs a higher angular momentum difference as evident from 
equation (|12|l . 

Via the angular momentum conservation the flows are 
more tightly constrained than in the original SB09 model. 
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They now depend on physical properties that are prin- 
cipally observable and connected to physical processes. 
But even without any observational constraints, the chosen 
parametrisation leaves less freedom, for example it is not 
possible to provide the total gas mass by radial flows alone, 
whereas pure infall is still an option. 

This treatment is extremely simplified and does in no 
way describe the whole complexity of the processes involved, 
ft would be naive too believe that the accreted material 
has some smooth and fixed angular momentum. We should 
hence bear in mind that the laid out formalism rather gives 
effective profiles resulting from all processes that may af- 
fect the gas before it mixes with the disc gas. Despite these 
limitations the result still provides useful information on 
these processes. We stress that the proposed coupling be- 
tween onfall and radial flows together with the constraint 
of an exponential density profile now completely determines 
the spatial structure of onfall and radial flows. As the ra- 
dial density profiles are kept fixed the mass distribution in 
each ring is not changed by the flows, but rather only by 
the continued accretion onto the galaxy. Hence, the radial 
flows only affect the distribution of the produced metals. 
Instead of using a preset formula for the flows throughout 
the disc these are derived consistently from the requirement 
of angular momentum conservation and the total required 
mass. 
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4 MODEL BEHAVIOUR 

This part presents the results of the new prescription for cou- 
pling flows and onfall by the conservation of angular momen- 
tum. We start by reproducing the flow profile of the SB09 
model to show that the mechanism works in principle. Then 
we discuss the basic dependency of the main observable con- 
straint, the metallicity gradient, on the model parameters by 
showing illustrative examples. 

In Fig[T]the SB09 model is compared to the new scheme 
with parameters a = 1/2 , b = 0. All radial profiles are 
shown for the present day disc age of 11.7 Gyrs. We see an 
acceptable agreement for this simple choice of parameters. 
This demonstrates on the one hand that the flows described 
by the naive formalism of SB09 were not unphysical and 
on the other hand we have an excellent starting point for 
our analysis by accurately reproducing the results of the 
older approach. A minor discrepancy resides in the inner- 
most rings. It can be traced back to the huge relative an- 
gular momentum changes required for driving gas inwards 
in the central regions inhibiting proper radial flows, while 
the older approach simply set some fixed flow. Anyway this 
region is not of interest here due to a general lack of obser- 
vations and the inadequacy of any model without a reliable 
description of the Galactic bar. 

We proceed to further study the behaviour of the flows 
for different parameter values. In Fig. [2] the resulting pro- 
files for different choices of parameter a and b are shown. The 
average angular velocity of onfalling material v = 0.5 a + b 
will turn out to be an important quantity for the flows with 
direct consequences for the metallicity gradient of the simu- 
lations. For the models shown it is v = 0.7, 0.3, 0.6 and 0.6, 
and we clearly see how the flows decrease with v increasing. 

Thus, adjusting the parameter b mainly sets the overall 
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Figure 1. Upper panel: the rate of flow of gas over the annulus 
with radius R, for the SB09-model and the new inflow model with 
a = 1/2, b = 0. Lower panel: the rates of accretion from the IGM 
per unit area of the disc. 



levels of inflow across the Galaxy as shown in the second (red 
diamonds) and third (green squares) model and it affects the 
whole disc in a similar fashion. The effects of parameter a are 
strongest in the outer regions of the Galaxy as a determines 
the angular momentum of the outer disc accretion and hence 
dominates the inflow behaviour there. Comparing the second 
(red diamonds) and fourth model (blue points) we see that 
increasing the parameter a leads to a stronger decline of 
the flows towards the outer regions of the Galaxy whereas 
the overall level is less affected than by changing b. In this 
way the radial extent of a high flow regime is controlled by 
setting the parameter a. Finally it is interesting to consider 
the effects on the flows when varying a and b while keeping 
v fixed as for the third (green squares) and fourth (blue 
points) models. A higher a in this case implies lower angular 
momentum in the inner regions which leads to a stronger 
buildup of the flows, and a higher maximal flow velocity. In 
the outer regions the angular momentum is actually higher, 
and the flows consequently decline. At fixed average angular 
velocity increasing a increases radial flows in the centre of 
the Galaxy and decreases them in the outer regions. 

In the lower part of Fig. [2] we see the corresponding ac- 
cretion densities. Albeit flattening at lower onflow angular 
momenta, which support the inner disc by radial gas flows, 
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Figure 2. Upper panel: the rate of flow of gas over the annulus 
with radius R for new inflow model for the best fit model with 
a = 0.33 for b = 0.53, a = 0.2 for b = 0.2, b = 0.5 and b = 0.7 
and a = 0.2, a = 0.5, a = 0.8 for b = 0.2. Lower panel: the rates 
of accretion from the IGM per unit area of the disc. 



the accretion density is nearly exponential over most of the 
disc in all cases preserving the shape of the exponential sur- 
face density throughout the simulation. 

As mentioned before, it should be checked whether the 
flow velocities stay within reasonable limits: we show this 
in the upper panel of Fig. [3] for different parameter values. 
During the buildup of the Galactic disc, little gas in the disc 
faces a fulminant onflow leading to high radial velocities as 
demonstrated by the streak of high velocities at the earli- 
est times in the lower panel. During that period the flows 
are in some cases limited by the total mass content of the 
next outer ring. This issue is not overly important, since the 
behaviour at those early times does almost not affect the 
present day metallicity gradient. For parts of the parame- 
ter space (at low specific angular momentum) flow velocities 
can even get so high that they cannot be sustained at later 
times. Fortunately this is limited to otherwise uninterest- 
ing cases, where the onfalling gas has extremely low angular 
momentum. In that case this parametrisation fails and simu- 
lations with such parameter values should be discarded. For 
the parameter range of interest this behaviour is roughly 
limited to the first Gyr of the simulation as can be seen 
in the lower panel of Fig. [3] Anyway the model currently 
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Figure 3. Upper panel: Resulting flow velocities for models with 
a = 0.5 and b = 0.0, 0.3, 0.5 at times 3, 6, 9 and 12 Gyrs. Lower 
Panel: flow velocity map for the best fit model with a = 0.33 and 
b = 0.53 across the disc and the simulation time. 



does not provide a satisfactory description of early Galaxy 
formation. 

The high flow velocities seen during the early times in 
the simulation should not be confused with an inside-out 
formation of the disc. Flow velocities increase outwards due 
to the build up of flows throughout the disc and the masses 
in each ring decreasing according to the exponential density 
profile. 

As our goal is to constrain the flows by fitting the Galac- 
tic radial abundance gradient we need to understand how the 
gradient depends on the flows. In Fig.|4]the resulting [Fe/H]- 
gradient with linear fits in the region 4.0 kpc < R < 14.6 kpc 
is shown. As expected, an increase in onfall angular mo- 
mentum, which lowers the flows, also lowers the metallicity 
gradient. In addition we clearly see the effect of depletion 
in the outer regions of the disc. If higher flows build up 
throughout the disc, the outer regions grab more gas from 
the IGM, which leads to a general decrease in metallicity 
and a stronger peak in the central regions, where the met- 
als from stellar yields get advected. At the outermost ring 
the assumed metallicity of the fresh material dominates the 
local metallicity, since the model does not take into account 
any star formation beyond that point, so that the curves 
meet there. The second (red diamonds) and third (green 
squares) models again illustrate the effect of parameter b on 
the gradient. An overall lower flow level lowers the metallic- 
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ity gradient. Moreover, lower overall flows limit the penetra- 
tion depth of the depletion zone into the galactic disc which 
then leads to a sharper drop off in the outermost parts. The 
comparison between the second model (red diamonds) and 
the fourth model (blue points) shows the effects of chang- 
ing parameter a. A higher value a decreases the flows in 
the outer disc relative to the inner disc and hence causes 
the shown upwards bending of the metallicities. Comparing 
model 3 (green squares) and 4 (blue points) which have the 
same average angular velocity, we observe that the gradients 
agree fairly well, at —0.0799 dex/kpc and —0.0795 dex/kpc 
respectively, just the curvature is different. Consistent with 
the flow patterns discussed above, at fixed average angular 
velocity the model with higher a has a higher metallicity 
and a steeper gradient in the very centre which flattens out 
towards the solar radius and the outer parts of the Galaxy. 

Summarising the main points of this part we see that 
the model only works for sufficiently high angular velocity of 
the onfalling material corresponding to sufficiently low radial 
flows. However, flows exceeding our threshold of 8.2km/s 
produce unphysically large gradients and the metal advec- 
tion to the central regions would require significantly larger 
mass loss of the central regions than currently assumed. The 
threshold is hence of no interest for the parameter region un- 
der study. The abundance gradient behaves as one naively 
expects with some complications due to the outer boundary. 
Fortunately, these boundary effects do not seriously affect 
the regions of the disc with R < 17kpc. 



5 FITTING THE MODEL 

In this section we will fit the model to data to the Galac- 
tic abundance gradient. We start by comparing the simu- 
lations with the compl ete catalogue of Cepheids given by 
iLuck fc Lambert! (|20 1 if ) and references therein. They anal- 
ysed a total of over 400 Cepheids and show that the 398 
stars passing their quality cuts deliver a linear relationship 
of metallicities with Galactocentric radius. [Fe/H](i?) is best 
fit by [Fe/H] = (-0.062 ± 0.002)i?/ kpc + (0.605 ± 0.021). In 
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Figur e 5. [Fe/H] from Cepheid Data as given in lLuck &: Lambertl 
a201lh with both best fit models, the first obtained without an 
metallicity offset, the second with an offset by +0.14 dex. 

Fig.[S]we see the observed data points and the best-fit mod- 
els, the first fitted without and the second with a metallicity 
offset of +0.14 dex discussed below. 

The parameter space of the new flow model is explored 
on a grid consisting of 992 simulations equispaced on the 
(a + 6)-&-parameter plane. When interpreting the following 
plots it is useful to keep in mind that the case of constant 
angular velocity (a = 0) corresponds to the diagonal from 
the lower left to the upper right corner. To the left and over 
this line a < while to the right and below it a > 0. Lines of 
constant average angular velocity v = 0.5 a + b correspond 
to parallels to the diagonal from the upper left to the lower 
right corner. 

To assess the agreement between the data set and the 
simulations we computed 
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where yi, S im is the [Fe/H] value of the simulation at radius 

Xi. 

In Fig. |S] we plot xved 2 -values for the entire grid and 
parameter slices intersecting at the absolute minimum are 
displayed on the rim. Contours are shown on the lower panel. 
By looking at the contour plot the unique minimum at a = 
—0.233, b = 0.8666 can be identified. This corresponds to an 
average angular velocity of the onfalling material of v/V c = 
0.75 or v = 165km/s assuming V c = 220km/s. 

With this naive approach the fit results depend strongly 
on the detailed assumptions about chemical evolution, be- 
cause a change in the total metallicity will shift the mini- 
mum and alter its depth in our parameter space. The abso- 
lute metallicity is e.g. impacted by the specific stellar yield 
dataset in use or by a different IMF, which changes the ra- 
tio between lock-up and yields, or by the overall loss rate in 
Galactic winds, factors which are quite weakly determined. 
Aditionally the determination of stellar abundances may suf- 
fer from systematic errors of the order of 0.1 dex. Especially 
the Salpeter IMF has the largest mass of the known IMFs 
in use, which leads to a probable lock-up of material in low- 
mass stars, which in turn lowers absolute metallicity levels. 
To the contrary, relative abundances are known to high ac- 
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curacy and the gradient in the simulations is far more robust 
as well. 

So for this purpose it is far more reliable to treat the 
average metallicity as a free parameter. This still gives a 
unique best model and improves the fit for parameters away 
from this minimum. The quality of the best fit is compara- 
ble for both procedures as they differ mostly in the inner 
parts of the Galaxy where no data constrains the models. 
The best fit with an allowed offset shows a linear gradient 
across almost the whole Galaxy, whereas the model fitted 
without the offset flattens towards the centre. Effectively, 
allowing for an offset fits the form of the gradient resulting 
in linear abundance gradients over the whole disc for the 
best fit model. 

Fig. [7] shows the corresponding results when fitting the 
average metallicity of the simulations to the data. The off- 
sets used for the best fits are shown in Fig. [8] Generally 
higher flow velocities decrease the average metallicity in the 
simulations, but the offset is relatively small in the region 
where the best fits are achieved. The absolute minimum is 
now at a = 0.3333, b — 0.5333 with an offset of 0.14 dex, 



the simulations having the lower average metallicity. The 
average angular velocity is v/V c = 0.7 or v = 154km/s. 

Finally we note that the metallicity gradient of the sim- 
ulations is tightly correlated with the average angular veloc- 
ity v of the onfalling gas as seen in Pig. El This relation grows 
worse when using a different value, e.g. the flow velocity at 
the solar radius, instead of the average. Thus, the results 
lead us to consider the average angular velocity as an impor- 
tant quantity to understand the simulations. As discussed in 
the preceding section at fixed average angular momentum 
the models with higher a tend to have a slightly flattened 
metallicity gradient. This creates some scatter around the 
average relation depending on the precise values of the pa- 
rameters a and b which justifies fitting the models on the 
whole grid instead of only working with the average veloc- 
ities. But the qualitative behaviour of the simulations can 
clearly be extracted from this plot alone and the best fit can 
be determined up to an uncertainty of Av/V c = 0.07. 

In summary, the average velocity can be constrained to 
be 0.70 V c = 154km/s v s: 165km/s = 0.75 K . The 
average velocity would be a plausible value for gas from 
the corona. The corona presumably rotates slower than the 
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steepen at the same average angular velocity v of the on- 
falling material. Hence, the best fit should shift to higher 
values of v reducing the flows caused by onfalling material. 
However, the nature of these flows is drastically different 
from those caused by the onfalling material, which is not 
constant, but varies radially due to the coupling to the in- 
fall. 

Specifically, we add an additional constant inflow with a 
speed of v = 0.2, 0.4 and 0.6km/s to the simulations. Fig llOl 
shows the resulting metallicity gradients for the simulations. 
We see that the additional flow has a significant effect over 
the whole parameter range. For v ^ 0.4 the simulations with 
the additional flows show a considerably steeper gradient, ly- 
ing below the reference values (red squares). Below v = 0.4 
the simulations show a shallower gradient being above the 
reference values. However, this region is not important as 
the simulations there show an unreasonably high gradient 
that is not observed. Moreover, the flows in this regime are 
so strong that they are actually limited by the maximally 
allowed flows of the simulation. For these parameter values 
the flows saturate in the centre of the galaxy concentrating 
most of the produced metals in the innermost ring and the 
strong inflow across the whole disc transports low metallic- 
ity gas from the outskirts of the Galaxy towards the inner 
regions reducing the average metallicity and gradient there. 
Thus, we disregard the models in this parameter regime. 

In the region of interest, the differences roughly amount 
to Au/Vc = +0.07 for the simulation with an additional 
constant flow speed of 0.2 km/s, and the differences are cor- 
respondingly higher for the other simulations at Av/V c = 
+0.11 and Av/V c = +0.15 for 0.4 km/s and 0.6 km/s re- 
spectively. 

We conclude that additional radial flows across the 
whole Galactic Disc strongly influence the abundance gra- 
dients. But the flows due to onfall remain important for the 
gradient even in this case. The models considered here are 
extreme cases in the sense that they assume an additional 
constant inflow velocity that extends across the whole disc 
and consequently the effects on the gradient are relatively 
strong. Moreover, only the first model is actually within the 
estimates on the order of the flows caused by the considered 
effects, and the second and third model already have higher 
than expected flow levels. 



disc gas by about vi ag ^ cor « 75 km/s, which corresponds to 
v cor « 145 km/s when neglecting any vertical velo city gra- 
dients |Marasco et al.ll201ll ; iMarinacci et alj|201ll ). Appar- 
ently, this result depends on the gradient to be fitted and 
we cannot claim that this result would be unique as there 
are different ways of producing a gradient. Still, there seems 
to be a strong tendency to recover such values within this 
framework. 

5.1 Additional inflow 

As discussed before, both (turbulent) friction in the gas lay- 
ers and interaction of the gas with spirals and bars can drive 
a radial inflow which could be of the order of 0.3 km/s. To 
ascertain how strong this affects our results, we computed 
simulations in which a constant radial inflow velocity was 
added to the inflow model. Because this simply increases 
the overall flow levels, we generally expect the gradient to 



5.2 Effects of a different infall history 

The prescription for the infall used in the original model 
might also influence the results obtained from the flow 
model. In particular, the evolution of the infall determines 
the available amount of gas and the SFR history. With a 
non-linear SFR law increasing the mass directly creates a 
gradient due to the stronger metall production in the centre 
of the Galaxy. As in our model the radially flowing mass is 
a fraction of the infalling mass, lower infall also reduces the 
flow levels across the whole disc yielding lower metallicity 
gradients. In addition a lower gas mass lowers the maximal 
gas flows, which may lead to differences from the expected 
behaviour for very high flow velocities, and increases the 
effects of flows on the gradient, an effect studied below. 

In the standard model infall increases the gas mass in 
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Figure 10. Metallicity gradient of the simulations as a function of 
the average angular velocity of the onfalling gas, for the reference 
simulation and with an additional inflow with constant speed of 
v = 0.2, 0.4 and 0.6km/s for models 1,2 and 3 respectively 
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Figure 11. Metallicity gradient of the simulations as a function 
of the average angular velocity of the onfalling gas for the simu- 
lations with a fast or slow decline of the infall and the standard 
model 



the disc according to 
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M — — — e ' 1 + -t — e 

01 02 

with the parameters bi = 0.3 Gyrs and b 2 = 14 Gyrs. To test 
the dependence on the assumed infall history the parameter 
62 of the second exponential was changed to i>2 = 7 Gyrs and 
&2 = 28 Gyrs corresponding to a faster or slower decline of in- 
fall respectively. Differences between the simulations should 
be higher for the case of the faster decline because the orig- 
inal simulation already has a relatively slow decline and the 
exponential decay law. 

In Fig. [TT] the results of the simulations are shown as 
a function of the average angular velocity. The simulations 
with the faster decline shows a considerably lower metallic- 
ity gradient for almost all parameter values in agreement 
with their lower flow levels, lying above the standard model 
in the plot, with a maximal difference up to 0.02dex/kpc. 
For decreasing v, thus, higher flows, the difference to the 
standard model gets smaller. In this parameter regime, the 
standard model's flows start to be limited by the total avail- 
able gas mass, while this is not yet the case for the generally 
lower flows of this model. 

The simulations with a slower decline show metallicity 
gradients comparable to the standard model, but with a 
slightly higher gradient for most of the parameter values 
corresponding to generally higher flows. For the very highest 
flow velocities, the gradient is actually lower. However, this 
difference is quite small and may be explained by saturation 
of the flows in the inner region of the galaxy and higher flows 
in the outer parts reducing the average metallicity. 

Regarding the best fit to the observed data, these dif- 
ferences would translate at most into Av/V c = 0.1 in the 
region of interest. 

A more detailed analysis can again be performed on a 
grid of the parameter values, in the case of the fast decline 
model the best fit values are a — 0.3, b — 0.5 corresponding 
to v/V c = 0.65, for the slow decline we find a = 0.366, 
b — 0.533 corresponding to v/V c = 0.716, both quite close 
to the standard model values of a = 0.333, b = 0.533 and 
v/V c = 0.7. 



In summary doubling the decay time has a small effect 
of Av/Vc = +0.016, whereas halving it leads to a difference 
of Av/Vc = —0.05 compared to the standard model. 

5.3 Effects of different star formation rates 

A different star formation efficiency has an effect on the flow 
model by changing the amount of gas available in the disc 
and on the metallicity by changing the SFR history. The 
flow model is affected as changing the gas mass in the disc 
changes the fraction of onfalling to disc gas and the total 
angular momentum in the disc. For higher gas masses in the 
disc the effect of flows are reduced as the same amount of 
metals transported radially by flows is diluted in a higher 
amount of disc gas. Additionally, the simulations will show 
different behaviour in the high flow regime when the flows 
are limited by the total mass in the rings. 

In the standard model stars form according to 

d^star , yil.4 

Changing the efficiency of star formation (k) also changes 
the stellar mass in the disc and thus the metals produced 
in stars. To avoid this problem and to isolate the effect on 
the flow model we kept the star formation history of the 
simulations constant by also using different values for the 
masses Mi and M2 in the infall law. When increasing the 
star formation efficiency, the available gas mass was lowered 
and vice versa. Thus, mainly the gas mass in the disc is 
affected. However, the same amount of metals produced is 
then diluted in more or less gas, also changing the gradient. 

The standard model with k = 1.2 X 10~ 4 , Mi = 0.45 x 
10 10 M©, M 2 = 2.9 x 10 10 M© was compared to simulations 
with a higher efficiency and lower gas mass k — 3.167x 10 -4 , 
Mi = 0.11 x 10 10 M©, M 2 = 3.07 x 10 10 M© as well as 
simulations with a lower efficiency and higher gas mass k — 
0.3789 x 10" 4 , Mi = 1.5 x 10 10 M©, M% = 2.50 x 10 10 M©. 

We expect the simulations with the higher gas masses 
to have a smaller metallicity gradient because higher gas 
masses should decrease the radial flow velocities of the gas 
at fixed v. This is generally confirmed in Fig. [12] which shows 
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Figure 12. Mctallicity gradient of the simulations as a function 
of the average angular velocity of the onfalling gas for simulations 
with different star formation efficiency 



the resulting gradients in the simulations. The simulations 
with the lower efficiency and higher gas mass show a shal- 
lower gradient over most of the parameter values. But for 
very low flow velocities, the gradient is steeper than in the 
standard model. Here the effects of a higher gas mass dom- 
inate over the effects of the flows and steepen the gradient. 
This is best seen for v — 1.0, which corresponds to no flows, 
for which the low efficiency high gas mass simulation shows 
the strongest gradient, the standard simulation a somewhat 
smaller gradient and the high efficiency low gas mass simu- 
lation a nearly flat gradient. 

For the high efficiency simulation three regimes of pa- 
rameter space in which the simulations show different be- 
haviour compared to the standard model can be identified. 
The discussed low flow regime, where the effects of gas mass 
and dilution dominate. The intermediate regime, where the 
direct effects on the flow model become visible, and the sim- 
ulation has a steeper gradient corresponding to it's lower gas 
mass. And the very high flow regime, in which the gradient 
flattens slightly because the flows saturate in the centre of 
the galaxy and transport more low metallicity gas into the 
inner regions lowering the average metallicities. 

Again the differences in the region of interest are rela- 
tively small and give roughly Av/V c = 0.1, the lower effi- 
ciency simulation yielding a lower average velocity and the 
higher efficiency simulation giving higher values. 

We conclude that different assumptions for the star for- 
mation rate do influence our results moderately by changing 
the gas mass in the disc. Fortunately, the effects are rela- 
tively small in the region of the observed metallicity gradi- 
ent giving an uncertainty of Av/V c = 0.1 for the simulations 
presented here. 



6 THE ACCRETION PROFILE 

While the angular m omentum is in formidable agreement, 
M arasco et alj (|201lh also predict a specific accretion pro- 
file for their theory of supernova-driven gas accretion which 
greatly differs from that seen in our models. 

To test whether their accretion profile allows for a plau- 



sible chemical evolution, we computed simulations based on 
the SB09 model. Assuming the accretion profile, we nor- 
malise the total accreted mass in each timestep to be the 
need of the whole disc and distribute the infall over the an- 
nuli according to the profile. The remaining missing mass in 
each annulus is then provided for by radial flows. 

Results of this model are shown in Fig. [T3] As the peak 
of the accretion rate is farther out than the peak of the 
star formation rate, the inner regions of the disc need to 
sustain their mass by high radial inflows. These high flows 
cause a steep gradient to build up in the inner half of the 
disc (R < 9kpc). As the infall again falls off beyond 9kpc 
more of the need in each ring is provided by radial inflows, 
but without the diluting effect of the low-metallicity infall 
from the IGM the metallicity increases between 10 kpc and 
15 kpc. Beyond that radius both the low star formation rate 
and the inflow from the outermost rings serve to decrease 
the metallicity in the gas. 

There are two main shortcomings of the resulting abun- 
dance gradient. First of all the very steep abundance gradi- 
ent of —0.164 dex/kpc in the region R < 9 kpc. This is an 
immediate consequence of the high radial flows necessary to 
redistribute the infalling gas inwards to the maximum of the 
star formation rate in the Galaxy. If such an accretion profile 
without additional infall from a different source is adopted 
this effect is inevitable and strongly disag rees with the gra- 
dient from observational data. Eventhough lLuck fc Lambert! 
l|201ll ) find a somewhat steeper gradient of —0.103 dex/kpc 
in the inner disc R < 6.6 kpc, this still differs greatly from 
the model results. Secondly the increase of the metallicity 
between 10 kpc and 15 kpc is definitely ruled out by the 
data. However, it might be easier to remove this problem 
by adding a different source of infall that provides enough 
gas in the outer rings to reduce the metallicity but does 
not greatly affect the global accretion rate due to the small 
masses required. In the inner disc it would be difficult to 
solve this problem, because we are facing an actual mass 
deficit and not a mere gradient problem, so we cannot sim- 
ply blame it on t he bar. The most lik ely way out is that 
the assumption of iMarasco et al.1 (|201lh of a minimum alti- 
tude above the disc that a gas cloud has to reach in order 
to accrete coronal gas fails or that multiple supernovae can 
assemble more gas there. 



7 NON-CIRCULAR ORBITS AND ANGULAR 
MOMENTUM TRANSPORT 

As mentioned in the introduction angular momentum might 
also be transported by stars moving on non-circular orbits. 
We assume that a star is born at some radial position with 
the corresponding angular momentum of the gas. During 
its lifetime it may leave this radial region and eject mate- 
rial at a different radial position. Depending on whether it 
spends more time outside or inside its original radius this 
will cause inflows or outflows. Two processes that are in- 
cluded in the SB09 model are important for this process: 
churning and blurring. If one uses a simplified picture of 
stellar orbits showing radial oscillations around their guid- 
ing centre radius, churning changes the guiding centre ra- 
dius and the effects of the radial oscillation are described by 
blurring. Churning exchanges stars and gas between differ- 
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Figure 13. Upper panel: the rate of flow of gas over the annu- 
lus with radius R for the SB09-model with the accretion profile. 
Middle panel: the rates of accretion from the IGM per unit area 
of the disc. Lower Panel: Abundance gradient assuming the ac- 
cretion profile of Marasco et al. (2011) 



their guiding centre radius. If a star ejects material at such a 
point, its yields have a different specific angular momentum 
than the surrounding gas, which will cause in- or outflows. 
Generally a larger dispersion makes populations expand out- 
wards. This process depends on the shape of the potential 
and the radial density and velocity dispersion profiles, but 
in general there is a positive asymmetric drift that drives an 
inflow. 

The model with parameters a = 0.5,6 = 0is used to 
estimate this effect because this choice closely resembles the 
original SB09 model. At this point the resulting flows are 
only calculated and not actually implemented. This should 
not have any impact on our conclusions in this part as the 
mass distribution is held fixed by the onfall mechanism. For 
the same reason the model parameters a and b have no effect 
on the velocity profiles and are only given for completeness. 
If the flow velocities are comparable to those resulting from 
the onfalling gas, it would certainly change the chemical 
composition of the disc and need to be considered. 

In Fig. [14] the results are plotted in the form of a ve- 
locity map covering the complete simulation. In the inter- 
mediate disc regions the inflow (positive values) driven is 
strongest, since there is still significant stellar heat, the cir- 
cular geometry is less important and so the effective asym- 
metric drift takes its maximum toll. The small outflow seen 
in the inner rings should not be taken too seriously as it is 
likely caused by a numerical instability (the relative angu- 
lar momentum/radius change to reach the adjacent ring is 
far smaller for the outer neighbour, which is a mere effect 
of the grid size). This region is anyway uninteresting, since 
the bar is not modelled. The outer boundary produces a 
similar effect in the outermost rings where stronger inflows 
develop due to the fact that in the simulation no stars can 
come from beyond 20 kpc. The resulting flow velocities are 
of the order of v ~ 0.01 km/s, negligible compared to the 
flow velocities onfalling material causes, which were of the 
order v ~ 1.0 km/s. Therefore, their impact on the chemi- 
cal evolution will not be significant and the conclusions of 
the last section remain unchanged apart from slight correc- 
tions, which are certainly smaller than the expected errors 
and uncertainties due to other processes. This statement 
may change if stars can exhibit a significant friction with 
the surrounding gas. We know that the Heliosphere extends 
far beyond the planetary system of our Sun and that many 
stars produce visible bow shocks in the interstellar medium. 
In case this interaction with the surrounding medium can 
slow down the relative motion of the stars and is not only an 
expression of the momentum balance with the stellar wind, 
the calculated terms could increase significantly. 



ent annuli by means of a resonant interaction with transient 
spiral arms in the Galaxy. In this interaction total angu- 
lar momentum is conserved, the star actually lowers (when 
moving inwards) or increases (when moving outwards) its 
angular momentum. Flows associated with the spiral pat- 
tern are hence driven by the different radial profile of gas 
and stars, which are already in the model, and the gas los- 
ing a bit of angular momentum in shocks in the spiral arms, 
which is generally a minor effect. 

Blurring affects stars that do not move on circular or- 
bits and, hence, spend time at galactocentric radii away from 



8 INSIDE-OUT FORMATION 

In the following section we will implement a simple model of 
inside-out formation in the simulations. We will still assume 
an exponential gas density profile, but make the gas scale 
length a function of time. Then we have 

E g (R) = E g ,o(t)e- R/R ^ . (13) 

Accordingly, in each timestep of the simulation a new den- 
sity profile is calculated, from which the required amount of 
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Figure 14. Flow velocities due to gas ejected from dying stars 



inflow and onfall follows by the flow prescriptions 
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T,o(t) is then determined by the total gas mass in the disc 
at this time. If M gaa (t) is this total gas mass, then 



S (t) 



Mgas{t) 



where the sum runs over the annuli of the simulation. 

The evolution of the scale length itself is neither tightly 
constrained by observations nor by theoretical predications. 
Again for simplicity we will assume it to be linear, i.e. 



R g (t) = Ri + (R f - Ri)t/t f 



(14) 



where Ri, Rf and tf are the initial and final scale length 
and final time respectively. We will consider two cases, the 
first in which the gas scale length grows linearly over the 
whole duration of the simulations, the second in which the 
growth is limited to a time t / smaller than the present day 
disc age of 11.7 Gyrs and constant thereafter. Note that in 
the second case the stellar scale length will still continue 
to grow after the point when the gas scale length remains 
constant because the stellar population needs time to adjust 
to the new scale. 

In the reference model the gas scale length is constant 
in time at R g = 3.5 kpc corresponding to a stellar disc scale 



length of Rstar = 2.5 kpc = R g /1A taking the exponent 
1.4 in the star formation rate into account. T his value is 
in agr eement with recent studies of the MW bv lJuric et alj 
( 2008) who find (2.6 ± 0.52) kpc for the thin disc. As con- 
straints for the final and initial values we initialize the gas 
disc at Ri = 2.80 kpc and let it grow up to Rf — 3.70 kpc 
for tf = 11.7 Gyrs in the first and tf = 4.5 Gyrs in the sec- 
ond case. These values correspond to an initial stellar disc 
scale length of Rstar = 2.00 kpc, consistent w ith the value 
of (1.96 ±0.05) kpc found bv lBovv et ail (|201ll ) for old stars 
in the MW. The slightly higher final gas scale length is nec- 
essary when aiming to reach the stellar disc scale length of 
the reference model despite the lag of the stellar popula- 
tion behind the gas disc. Values significantly higher than 
the final sc ale length for the gas co nsidered here seem to 
be unlikely, iKalberla fc Pedes! (|2008l ) find a scale length of 
3.75 kpc for the MW. Moreover, we will see that the lag 
introduced by former star formation is so strong that a pro- 
longed episode of significant gas disc growth seems hardly 
reconcilable with the currently accepted observations. The 
formation time for the second model was chosen such as to 
roughly reproduce the stellar scale length of the reference 
model at 11.7 Gyrs. Then the stellar disc will have grown 
by 25% over the simulation time, as found for disc galaxies 
since 2 = 1 (|Truiillo fc Pohlenl 120051 ; iMunoz-Mateos et all 

Eon]). 

We consider these models to investigate the general ef- 
fect of an inside-out formation on the simulations and our 
results for the angular velocity of the onfalling material. The 
treatment is certainly oversimplified and is not meant to 
model the details of the actual formation processes or give 
the correct time dependence. However, even with such a sim- 
ple toy model we will at least see the most important trends. 

As the inside-out formation changes the star formation 
history, to make the simulations better comparable to the 
reference model, we adjust the mass parameters Mi and M 2 
in the onfall law given by 



»V M i - 
M = -r- e 

bi 



t/bj 



, M 2 _ 



t/b 2 



to yield similar results for the star formation at the solar 
radius. Recall that b\ was a short time scale to model the 
build up of the Galaxy and fc 2 a longer time scale respon- 
sible for sustained star formation. Specifically, the reference 
values are Mi = 0.45 x 10 10 M Q and M 2 = 2.9 x 10 10 M , 
the first inside-out simulation uses Mi = 0.75 x 10 10 Mg, 
M 2 = 2.85 x 10 10 M Q and the second M a = 0.70 x 10 10 M , 
M 2 = 2.75 x 10 10 M . The increased mass for the short 
timescale onfall serves to produce the early peak in the SFR 
and the mass for the long time scale is decreased to account 
for the additional gas mass due to a higher scale length at 
later times. 

The corresponding star formation rates as functions of 
time at R = 7.5 kpc are shown in Fig 1151 As explained the 
parameters are chosen such that they roughly show a compa- 
rable evolution. In the case of the first inside-out simulation 
this is not entirely possible due to the ongoing growth of 
the gas disc over the whole time. The lower gas scale length 
redistributes the mass towards the centre of the Galaxy and 
to keep the SFR high enough farther out would require con- 
siderably higher total masses in the models. We decided to 
keep the final SFR comparable which leads to a deficit at 
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Figure 15. SFR at R = 7.5 kpc over the simulation time for the 
reference model and the two inside-out scenarios 



Figure 16. Gas and stellar disc scale lengths over the simulation 
time for the reference model and the two inside-out scenarios 



earlier times as the scale length of the gas disc is smaller. 
For the second inside-out simulation the agreement is qual- 
itatively better. In both cases the final and peak values of 
the SFR agree fairly well, but the peak is shifted towards 
earlier times. Without changing the functional dependence 
of the infall law itself a better fit can hardly be achieved. 

Next we consider how the gas scale length affects the 
stellar scale length of the simulations shown in Fig 1161 The 
first model shows a linear increase of the stellar scale length 
from the initial value of 2.0 kpc to a final value of 2.30 kpc 
failing to reach the reference value. Continuing the simu- 
lations for another 1.8Gyrs giving a gas and stellar scale 
length of 2.34 kpc and 3.84 kpc respectively is still not suffi- 
cient to remedy this problem. Looking at the second model 
we see a growth of the stellar scale length from the initial 
2.0 kpc to the final 2.5 kpc. Here the growth is linear up to 
a short time after the gas disc stops growing at 4.5 Gyrs, 
and then slows down but does not stop. In both simulations 
stars born at earlier times and smaller scale lengths survive 
to the present time and consequently reduce the scale length 
of the total stellar population. As the SFR peaks during the 
first 3 Gyrs stars born at that time contribute heavily at all 
later times. 

These problems can certainly be solved by assuming 
a different functional form for the time dependence of the 
gas scale length. However, this can only be accomplished 
by increasing the scale length at early times to reduce the 
contribution of stars born at low scale lenghts. The second 
model serves to illustrate how it then becomes possible to 
reach the reference value for the scale length at the final 
time. To achieve the same in the first model would require 
unrealistically high values for the final gas scale length of 
4.3 kpc in conflict with observations. 

Finally we turn to consider the effect of the inside-out 
formation on the metallicity gradient. Generally, an inside- 
out formation can steepen the gradient. This expectation is 
only partly true for our simulations as seen in Fig 1171 For 
the first simulation the gradient is indeed steeper for all the 
parameter values considered. As the gas has a shorter scale 
length over long parts of the simulation compared to the 
reference simulation, more stars and consequently metals 
are produced in the inner than in the outer parts of the 
Galaxy adding to the gradient produced by the flows. For 
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Figure 17. Metallicity gradient of the simulations as a function of 
the average angular velocity of the onfalling gas for the reference 
model and the two inside-out formation scenarios 



the second simulation two competing effects influence the 
metallicity gradient. With little to no flows (0.5a + 6 > 0.6) 
the inside-out simulation has a steeper gradient, in this case 
the initial formation period and the gradient created by a 
shorter scale length of the gas dominates over the effects of 
the flows. But in the parameter regime of high flows, the 
inside-out formation actually reduces the flows as the inner 
regions of the galaxy already have sufficient mass due to a 
shorter scale length and do not need high radial flows to 
sustain their mass profile. 

In summary, an inside-out formation scenario has no- 
table effects on the metallicity gradient, the first case would 
approximately yield Av/V c = 0.1. But the interplay between 
radial flows and the formation process complicate matters in 
such a way that the gradient does not always steepen as in 
the second simulations. In particular, we note that a simple 
linear growth of the gas scale length amplifies the effects of 
the inside-out formation and does not seem plausible consid- 
ering the high final gas scale length necessary. An inside-out 
formation that mainly occurs at early times and slows down 
has considerably smaller effects on the metallicity gradient. 
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9 ORIGIN OF INFALL 

So far we have discussed the implications of the coupling 
of radial inflows to an assumed accretion onto the Galactic 
disc and explored the behaviour of the model under dif- 
ferent assumptions on the parameters and underlying laws. 
But the question of the origin of the accreted material re- 
mained open. Different theories exist and our simulations 
show values for the kinematic properties of the onfalling ma- 
terial comparable to those predicted by an interaction of gas 
clouds with a rotating corona. However, it is widely believed 
that extended gaseous ou ter disc like that of the Milky Way 
|Kalberla fc Pe des 2008), are the dominant reservoir and 
source to nourish the inner disc's star formation (see e.g. 
ISancisi et al.|[2008l ') . So far it is not clear how this gas should 
loose its angular momentum to move inwards, but assum- 
ing such a mechanism was found, the consequence would be 
rather high radial flows of gas through the edge of the stellar 
disc. The implications of these radial flows on the chemical 
evolution of the Galaxy have not been considered before, 
but might be used to put limits on the possible accretion 
rates. 

As we saw before, by the intimate relationship between 
radial flows and radial abundance gradients shown e.g. in 
Fig. [9] the observed abundance gradient imposes strong lim- 
its on the maximum flows and hence on the maximum 
"through-disc" accretion of our Galaxy. 

To estimate the effect of inflow from such extended gas 
layers we compare the total accreted masses during the past 
3Gyrs, i.e. in the time from 8.7Gyrs up to 11.7 Gyrs by 
onfall onto the disc and by inflow into the outermost ring. 
We take the total mass accretion in the outermost ring as a 
good measure for the through-disc inflow. 

The sum of accreted mass and inflowing mass is in- 
dependent of the flow parameters as in our implementation 
flows do not change the mass profile of the disc. But as higher 
flows build up throughout the disc the relative importance 
of radial inflow into the edge of the disc should increase. In 
Fig.[TH]the mass fraction of inflowing mass to accreted mass 
for the last 4Gyr is shown for the reference simulation, the 
two inside-out scenarios from the preceding section, a model 
with a steeply declining infall rate, and the model with an 
additional fixed flow velocity of 0.6km/s. Recall that the 
first inside-out scenario had a continuously growing gas disc 
scale length, while the second only has a growth period in 
the first 4.5 Gyrs, and both models end up with a higher gas 
scale length than the reference model. The infall model has 
a faster exponential decay of the mass of infalling material, 
specifically a decay time of 7 Gyrs instead of 14 Gyrs. For all 
these models the mass fraction drops as v increases and flow 
levels across the whole disc and at the edge decrease. The 
inside-out simulations generally show a higher mass fraction 
than the reference simulation. As the gas scale length at the 
end of the simulation was in both cases higher than in the 
reference model, the mass in the outermost ring is higher 
as well and more mass flows in at the edge of the Galaxy. 
Additionally, the inside-out growth of the disc causes higher 
flows in the outer regions of the Galaxy and, thus, a higher 
flow into the last ring. The second inside-out scenario is 
again closer to the reference value as most of the growth oc- 
curred before the time-frame over which we add the masses, 
and only the effect of the higher gas scale length is rele- 
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Figure 18. Upper panel: the fraction of total inflow mass to 
accreted mass on a logarithmic scale as a function of average 
angualar velocity v, for the reference model, the two inside-out 
formation scenarios, the model with a fast decline of infall and 
the model with an additional fixed inflow. 



vant in this case. The model with a faster decline of the 
infall rate shows a smaller mass fraction as total flow levels 
are reduced compared to the reference simulation. Finally 
the model with an additional constant flow velocity shows 
a considerably higher mass fraction for most parameter val- 
ues. This is easily understood as this model generally has 
higher flow levels and always drives a radial flow in the out- 
ermost ring, which is not true for the models without this 
additional flow. Again the lowest values of v should not be 
considered for this model. 

To compare the simulations, note that the model with 
additional fixed inflow should be offset by roughly Av/V c = 
+0.15 and the model with a fast decline of the infall by 
Av/Vc = —0.05 to yield similar metallicity gradients. At 
v = 0.5 the first inside-out scenario yields a contribution 
of inflow of 11 % to the total accreted mass, for the other 
simulations the contribution at comparable values of v is 
even lower. This value corresponds to a metallicity gradi- 
ent shallower than —0.095 dex/kpc for the reference model. 
Even considering some of the possible uncertainties dis- 
cussed above v >= 0.5 seems to be a reasonable lower bound 
on the average angular velocity of onfalling material. 

We conclude from these results that radial accretion 
from the outer disc gas reservoir through the Galactic disc 
should be relatively small. Considering a moderate lower 
bound of v >= 0.5, the radial accretion contributes less 
than 11 % to the total accreted mass in the models consid- 
ered here. A comparable participation (fraction > 10%) of 
the outer gaseous disc in feeding the disc would create ra- 
dial abundance gradients far in excess of the observations. 
A mechanism for the material to lose most of its angular 
momentum and to distribute the material over the whole 
disc without causing flows within the disc, i.e. flowing above 
the disc and without participating in the advection of met- 
als towards the centre would be required to be in agreement 
with the observed metallicity gradient. Yet, currently such 
a construction appears artificial and lacks any observational 
evidence. 
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10 CONCLUSIONS 

While recent works concentrated mostly on describing the ef- 
fects of inflows of a certain preset velocity, we coupled the in- 
flow process to conservation of angular momentum between 
the accreted material and the gas already present in the 
disc. In this picture, the angular momentum of the disc gas 
is diluted by freshly accreted material with lower specific an- 
gular momentum driving a radial flow through the disc. This 
coupling is used to derive the spatial structure of the radial 
flows consistently within the model. The resulting flow and 
onfall pattern is not known from the beginning, but rather 
a prediction of the model. Due to its physical motivation 
the new prescription is preferable to the older prescription 
of SB09 and also to models with a preset inflow. 

We have shown that the chosen prescription can pro- 
duce a vast range of radial metallicity distributions in disc 
galaxies and studied the dependence of the observables on 
different assumptions about the qualities of the inflow. 

We fitted the model to observation s of the present 
day m etallicity gradient of our Galaxy bv lLuck fc Lambert! 
(2011). This comparison constrained the average rotational 
velocity of the onfalling material to be 0.7 ^ v/V c ^ 
0.75 corresponding to 154km/s ^ v on j a ii ^ 165km/s as- 
suming a constant V c = 220km/s. This value compares 
well with the predictions for an accretion of gas from the 
corona (when neglecti ng vertical velocity gradients) at about 
Vcoro na ~ 145km/s (|Marinacci et al.ll201ll ; [Marasco et al.l 
201lf). Howeve r , the specific accretion profile predicted by 
Marasco et al.l (|201 lh is not compatible with our simula- 
tions, since the rather sharp peak in Galactocentric radius 
would drive flows that create a radial metallicity distribu- 
tion at odds with the data, if not some other process fills 
the gap especially in the inner disc regions. 

One major uncertainty is how much angular momen- 
tum is given up by the gas itself through various channels. 
We studied a so far neglected process where the ejection 
of material from stars moving on non-circular orbits dilutes 
the angular momentum via their asymmetric drift. If there 
is not an unknown amplification, this momentum exchange 
has only a very small impact on radial flow velocities of the 
order of vji ow ~ 0.01 km/s and can be safely neglected. More 
of a concern are the well-known processes of gravitational in- 
teraction of the gas with spiral arms or the bar and possibly 
friction withi n the g as layers. Considering the estimates by 
lLacev fc~F all ( 1985) they might play a role in additionally 
changing the gradient. We have considered models in which 
a fixed constant flow velocity was added to our flow model to 
describe part of the uncertainty of these processes. We con- 
clude that flows at that level certainly change our results 
quantitatively if they extend across the whole Disc and per- 
sist for the whole age of the Galaxy, but for values within the 
estimates the total effects remain small of order 0.05 v/V c . 
However, the local structure and gas dynamics may well be 
affected considerably without affecting the global flow levels 
and the gradient. 

Resonances stable in galactocentric radius are another 
source of uncertainty to the flow pattern. Any significant 
influence, however, should show up as an excursion in the 
radial abundance profile and up to now is unobserved, mak- 
ing such a bias a minor concern. 

Radial migration (jSellwood fc Binnevlliool l has only a 



minor influence in flattening the inner disc abundance gra- 
dient slightly via mixing. Even larger amplitudes are easily 
compensated by the inflow. 

These results were shown not to depend strongly on 
the resolution or the churning amplitude. Additional depen- 
dence on model parameters, specifically the assumed infall 
history and the star formation efficiency, was found to be 
relatively small. 

Declining star formation rates tend to reduce the gradi- 
ent in our models. As the radial flows are directly coupled to 
the accreted material which determines the star formation 
rates, a declining SFR also means decreasing flows. For the 
specific model parameters, a very fast decline of the SFR 
shifts the gradient to lower values, but a slower decline has 
relatively little effect in steepening the gradient. 

The star formation efficiency mostly affects the simula- 
tion by controlling the total amount of gas in the disc. For 
different levels of radial flows, the effects on the gradient 
are qualitatively different, but for the region corresponding 
to the observed data a higher gas mass tends to flatten the 
gradient, hence needing stronger flows and a lower angular 
momentum of onfalling material to give the correct value. 

We also examined the impact of inside-out formation. 
Generally inside-out formation raises the ratio of metal pro- 
duction to infall of fresh material in the inner disc relative 
to the outer disc facilitating a steeper metallicity gradient. 
However, a simple linear growth of the gas scale length over 
the complete age of the Galaxy is strongly limited by obser- 
vations, and a model in which the gas disc mainly grows in 
the early times of the simulations is more adequate in order 
to reproduce the stellar disc scale length measured today. 
In this latter case the effects on the present-day metallic- 
ity gradient are greatly reduced. It was impossible for us to 
achieve the measured abundance gradient just by inside-out 
formation and so conclude that radial flows are essential for 
explaining the observations. 

Independent from these considerations we showed that 
the large amounts of HI gas in extended layers beyond the 
stellar disc cannot feed the majority of fresh material re- 
quired to the inner disc, since this would imply too large 
radial flows and a too steep abundance gradient. 
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and Vfiow <4.1km/s, 8.2km/s, 16.4 km/s for the high, nor- 
mal and low time resolution simulations. Tightening the ve- 
locity limit can artificially improve the fits in the low in- 
fall angular momentum region, where the flows are too fast, 
but in the interesting region has virtually no impact. We 
demonstrate this in Fig. and Fig. 1201 While both coarse 
resolution simulations show some minor deviation from the 
general picture, using higher resolution apparently brings no 
further benefits. 

To test the effects of the churning on the results simu- 
lations with a — 0, b — 0,0.1, ■ ■ ■ ,0.9, 1 were computed at 
identical resolution with and without churning. In its current 
implementation churning also involves the cold gas phase 
and directly mixes the gas apart from the stars and hence 
their yields being scattered over different radii. The results 
are shown in Fig. 1211 The only significant difference happens 
in the case (b — 1.0), where there are no systematic radial 
gas flows, which otherwise dominate any effect from churn- 
ing. In other cases, churning mixes especially the inner disc 
and by this mildly lowers the inner disc radial abundance 
gradi ent. If churning was hence significantly stronger than 
set in lSchonrich fc Binnevl l|2009l ). we need a bit lower inflow 
velocities or vice versa a stronger radial flow to obtain the 
same gradient. The effect is close to being negligible, how- 
ever. We also note that in the formulation we used, churning 
gets stronger for coarse spatial resolution and fine time res- 
olution. 



11 APPENDIX: RESOLUTION EFFECTS 

The resolution of the simulation has several impacts on the 
results, for example by "numerical mixing" on the flow, 
when the flow does not cross the width of one annulus 
in the simulation per timestep and hence the gas thought 
to be on a ring boundary mixes with the entire ring. For 
our tests we vary both the spatial and the time resolution 
separately by computing in addition to the standard pa- 
rameters AR = 0.25 kpc, At = 30Myr a coarse grid of 
models with AR = 0.1 kpc and AR = 0.5 kpc as well as 
At =15 Myr and At = 60 Myr. Further resolution effects 
result from the maximal flow velocity limit given by the 
grid size Vfi ow < AR/ At, which has some limited influence 
in the earliest timesteps. For the resolutions used this cor- 
responds to Vfi ow < 3.3 km/s, 8.2 km/s, 16.4 km/s for the 
higher space, normal and lower space resolution respectively 
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Figure 21. Resulting abundance gradients for models with different space resolutions and with or without churning, the flow parameters 
are a = and b is given at the upper left corner of the plots 



